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Heterostructures of spin-orbit coupled materials with s-wave superconductors are thought to be capable of 
supporting zero-energy Majorana bound states. Such excitations are known to obey non-Abelian statistics in 
two dimensions, and are thus relevant to topological quantum computation (TQC). In a one-dimensional system, 
Majorana states are localized to phase boundaries. In order to bypass the constraints of one-dimension, a wire 
network may be created, allowing the exchange of Majoranas by way of junctions in the network. Alicea et al. 
have proposed such a network as a platform for TQC, showing that the Majorana bound states obey non-Abelian 
exchange statistics even in quasi-one-dimensional systems. 1 Here we show that the particular realization of non- 
Abelian statistics produced in a Majorana wire network is highly dependent on the local properties of individual 
wire junctions. For a simply connected network, the possible realizations can be characterized by the chirality 
of individual junctions. There is in general no requirement for junction chiralities to remain consistent across a 
wire network. We show how the chiralities of different junctions may be compared experimentally and discuss 
the implications for TQC in Majorana wire networks. 

PACS numbers: 03.67.Lx, 03.65.Vf, 03.67.Pp, 05.30.Pr 



I. INTRODUCTION 

Quantum computation, based on the creative use of the fun- 
damental resources of quantum mechanics, promises expo- 
nential speed up of several classically intractable computa- 
tional problems. However, since quantum states are extremely 
susceptible to external perturbations, maintaining their coher- 
ence in the presence of environmental interactions is the fore- 
most challenge in any quantum computation architecture. In 
the recently proposed scheme of topological quantum com- 
putation (TQC), 2 3 the environmental decoherence problem is 
confronted by encoding quantum information in an intrinsi- 
cally non-local way, making it essentially immune to any lo- 
cal perturbation due to the environment. A growing class of 
theoretically predicted quantum many-body states, character- 
ized by excitations with non-Abelian statistics (non-Abelian 
anyons), allow such non-local encoding of quantum informa- 
tion. 

Exchange statistics^ is the description of how a many -body 
wave function transforms under the unitary transposition of 
any pair of quantum particles. The simplest examples of this 
transformation are those associated with bosonic (multiplica- 
tion by 1) and fermionic (multiplication by —1) statistics. In 
(2 + 1) -dimensions, where simple permutation of the coordi- 
nates and actual exchange of the quantum particles are not 
necessarily equivalent, the bosonic and fermionic statistics 
can be generalized to anyonic statistics. In Abelian anyonic 
statistics, the many body wave function can pick up any phase 
between and tt under pair-wise exchange of the particles, 
which are now called anyons. In (2 + 1) -dimensions, if the 
many-body ground state wave function happens to be a linear 
combination of states from a degenerate subspace, a pair-wise 
exchange (braiding) of the particles can also unitarily rotate 
the wave function in the ground state subspace. In this case, 
the braiding statistics can be non-Abelian 23 . The state rota- 
tion produced by braiding may be exploited to manipulate the 



quantum information stored in the ground state manifold, pro- 
ducing quantum gates that may be used for computation!^ 
Because of the non-local storage of information within the 
ground-state subspace, TQC using non-Abelian excitations is 
intrinsically fault-tolerant. This intrinsic fault tolerance at the 
hardware level holds considerable promise for the future suc- 
cess of quantum computation. 

Non-Abelian quantum systems in the so-called Ising topo- 
logical clasP are characterized by topological excitations 
called Majorana fermions. In some topological supercon- 
ducting (TS) systems, 9 Majorana fermions arise as non- 
degenerate zero-energy excitations bound to vortices of the 
superconducting order parameter. The second quantized 
operators, 7,, corresponding to the Majorana excitations 
are Hermitian, jj — 7^. This is very different from 
the ordinary fermionic (or bosonic) operators for which 
Ci cj. Therefore, each Majorana particle can be re- 
garded as its own anti -particle. 10 Majorana particles have 
been predicted to occur in some exotic many-body states 
such as the proposed Pfaffian states in the filling fraction 
v = 5/2 fractional quantum Hall (FQH) s ystem p^ spin- 
less chiral p-wave superconductors/superfluids,^ 2 ^ and non- 
centrosymmetric superconductors.^^ Recently, there have 
been proposals for manufacturing the necessary conditions for 
Majorana fermions by constructing heterostructure systems in 
which an s-wave superconductor is placed in proximity to the 
surface of a 3D strong topological insulator (Tlp^or a semi- 
conductor thin film with Rashba 17 18 or Dresselhaus^-type 
spin-orbit (SO) coupling. 

Following this, ID versions of these systems have been 
proposed! 2 ^" 2 ^ In ID it is not necessary to have vortex states 
for Majorana fermions to occur. Instead, Majoranas appear as 
zero energy modes trapped at the phase boundaries between 
topologically superconducting (TS) and non-topological su- 
perconducting (NTS) phases of the system. In effect, these 
would be practical realizations of a ID lattice model shown 
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earlier by Kitae\E! to contain such Majorana fermion end 
states. It is somewhat problematic for TQC applications that 
braiding operations are unavailable in one dimension. Al- 
though in principle the Majorana bound states may be moved 
by gating the system to allow expansion and contraction of 
the topological regions, 1 there is no room for the Majoranas 
to pass around one another in ID. This problem has been re- 
solved by Alicea et all 1 by the introduction of a so-called T- 
junction in a network of quantum wires, which allows one 
to exchange the end-state Majorana modes using a junction 
of three quantum wire segments. In this way, their proposal 
makes a step from a locally one-dimensional system to a glob- 
ally two-dimensional one, in which braiding is now possible. 

The non-Abelian statistics of Majorana fermions on a quan- 
tum wire network on a superconducting substrate is not obvi- 
ous given that the original arguments for non-Abelian statis- 
tics in p-wave superconductors relied on the Berry phase accu- 
mulated from taking a Majorana fermion around a vortexP^In 
their work, Alicea et. al have shown 1 how non-Abelian statis- 
tics arises in the wire network by approximately mapping the 
system to a ID lattice model similar to the one considered 
by KitaevP^l We reproduce this result in a more general set- 
ting, and show that the form of the braiding statistics actually 
implemented during an exchange of Majorana bound states 
within a wire network is dependent upon local characteristics 
of the wires and their junctions. 

We begin in Sec.|II]by formulating the problem of exchang- 
ing a pair of Majorana fermions in a fashion that is inde- 
pendent of the underlying wire network. First, we show that 
given a pair of Majorana operators the signs acquired by them 
during an adiabatic exchange must be unique and consistent 
with non-Abelian statistics. In other words, if 71 and 72 are 
the 2 Majorana fermions being exchanged, then the result of 
such an exchange is 71 —> A72, 72 — > — A71 and A 2 = 1. 
The result can be described in terms of a braid matrix writ- 
ten as U = e lA7172 . The braid matrix U associated with a 
given exchange operation will be shown to be uniquely deter- 
mined from the microscopic parameters of the quantum wire 
network. By reversing the trajectory of the exchange opera- 
tion the sign of A is also reversed. In the remainder of the 
paper, we elucidate the above in the context of a quantum 
wire network. We show in Sec. Ill that the phase acquired 



by the Majorana fermions when they are exchanged through a 
junction is determined by a local characteristic of the junction 
(the junction chirality) that is independent of the exact path 
taken by the Majoranas, as well as the locations of the TS and 
NTS regions. We discuss the implications of this character- 
istic (and particularly the possibility that it will be different 



for different junctions within the same network) in Sec. IV A 



before extending our analysis to multiply connected junctions 



in Sec. IV B There we see that in addition to the chirality 
of junctions there is another representation-invariant quantity 
relevant to the braiding of Majoranas, i.e. the phase acquired 
when a Majorana is transported around a loop in the network. 
This 'loop factor' completes the description of the transfor- 
mations produced by Majorana motion in wire networks. Our 
analysis appeals only to the general necessity for consistency 
in the effects of Majorana motion, as well as a few simple as- 



sumptions about the network set forth in Sec. Ill With these 
assumptions satisfied, the results apply to any network of ID 
wires supporting Majorana fermions at phase boundaries. 



II. MAJORANA FERMION EXCHANGE IN THE 
HEISENBERG REPRESENTATION 



The topologically degenerate subspace of states of a sys- 
tem of topological nanowire segments proximity-coupled to a 
superconductor may be manipulated via an adiabatically time- 
varying Hamiltonian. Such operations result in changes of ex- 
pectation values of the various observables composed of prod- 
ucts of Majorana fermion operators. These expectation values 
can be computed equivalently in both the Schrodinger and 
the Heisenberg representations. Therefore the non-Abelian 
statistics generated by braiding operations can be studied by 
analyzing the time-dependent Majorana operator jj(t) = 
U (f)ljU(t) in the Heisenberg representation. Here U{i) is 
the unitary time-evolution operator. 

In this section we will show that in a general supercon- 
ducting system an exchange of one pair of Majorana fermions 

71 and 72 in the Heisenberg representation is described by a 
non-Abelian braid matrix U(t) — e 4 *e ±7r7172 / 4 , where r is a 
time after the exchange is complete. Such a non-Abelian uni- 
tary transformation leaves all Majorana fermions other than 
71.2 unchanged. The transformation U interchanges 71 and 

72 with a relative — sign such that 71 (r) = 72(0) and 
72 M = -7i(0) or7i(r) = -72(0) and 72 (t) = 71 (0). Be- 
low we show how an exchange of Majorana fermion modes 
will lead to such a relative — sign in a general setting. 



A. Uniqueness and reversibility of exchange transformation 

We start by describing the process of exchanging a pair of 
Majorana fermions in terms of the underlying BCS Hamilto- 
nian in the Heisenberg representation. Consider first a Hamil- 
tonian Hbcs which has a pair of Majorana solutions 71 and 
72. Since 7^ are zero energy Majorana solutions they com- 
mute with the Hamiltonian ([HbcStIj] — 0) and they are 
self-adjoint (7J = 7^). In order to exchange the pair of Majo- 
ranas it is necessary to vary Hscs{t) adiabatically in time 
in a certain time interval [0,r]. The Hamiltonian is taken 
to be static before and after this interval. Thus, one can de- 
scribe the states before t < and after t > r by eigen- 
states of Hbcs(Q) and Hbcs{t). Moreover for exchange 
operations, we will require that the Hamiltonian at the end 
of the operation Hbcs(t) be the same as that at the be- 
ginning Hbcs{0) = Hbcs( t )- Since the evolution was 
adiabatic, zero energy Majorana operators at times t < 
evolve into zero energy operators at time t > r such that 
hj(t),H BCS (t)] = hj(r),H B cs(0)} = Ofort > r. If 
the time variation of the Hamiltonian is such that it physically 
exchanges the positions of the localized Majorana solutions in 
the time interval [0, t], then it follows that 



72(1") = si7i(0), 7i ( T ) = ^272(0), 



(1) 
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where s±.2 are constants. As a result of the application of the 
adiabatically time-varying Hamiltonian Hbcs(® < t < t), 
the Majorana operators will evolve according to the Heisen- 
berg equation of motion jj(t) = i\HBCs{t) iljit)]- The so- 
lution to this equation can be written formally in terms of a 
unitary operator U(t) as 



lj (t) = uHt) lj (o)u(t) 



(2) 



where U(t) is the time-ordered exponential U(t) = 
Te~ % $a dTH Bcs{r) _ p rom £q j^J it is clear that jj remains 
Majorana for the entirety of the time evolution (i.e. Jj(ty = 
7j(<)). Considering the square of Eq.|2]we find that 7j(t) 2 = 
1 since 7j(0) 2 = 1. Applying this relation to the Majo- 
rana transformation equation (Eq. 1 1, we find that s 2 = 1, 
so Sj = ±1. It follows from Eq. 2 that 71.2(1"), and there- 
fore S12, can be uniquely determined from the relevant BCS 
Hamiltonian. Using arguments analogous to the ones used 
above, one can show that if the trajectories are reversed i.e. 
by replacing Hscsit) by Hbcs{ t — t) the values of S\ and 
S2 are interchanged. Thus the exchange operation can be de- 
scribed in terms of a unique unitary operator U(t) which we 
will refer to as the braid matrix such that 



s l7l (0) = C/(r)t 72 (Q)[/(r) 
s 272 (0) = C/(r) t 7l (0)[/(T) 



(3) 



with the time-reversed braid matrix being described by 
U(t) = W{t). 



B. Non-Abelian statistics 

The above argument only shows that Sj = ±1. How- 
ever as discussed above, non-Abelian statistics implies a rel- 
ative — sign between the final Majorana fermions such that 
s%S2 = — 1. In this paragraph, we show by contradiction that 
this follows from the conservation of the fermion parity sym- 
metry obeyed by the BCS Hamiltonian. Suppose U = U(t) 
is the unitary operator that exchanges a pair of Majorana 
fermions 71 and 72. Suppose, for the sake of argument, that 
71 , 72 do not pick up a (relative) — sign under U. U then 
transforms the neutral fermion operator eft = 71 + 172 into 
IPdJU = ±id. Consider now the action of U on the ground 
state of Hbcs = -Hbcs(O) = Hbcs{t)- Since 7^ commute 
with the Hamiltonian Hbcs, so do S ,d and the number oper- 
ator d) d. Thus the ground state can be taken be a simultaneous 
eigenstate of d)d and Hbcs- Furthermore if Y$>) is a ground 
state so are d\^f) and d^ From here it is straightforward to 
see that the ground state has a 2-fold degeneracy: namely the 
empty state |0) and |1) = ef^O). Applying d) to U\0), where 
|0) is the empty state, we see that ct<U\0) = Ud\0) = 0, so 



U\0) = k\1) = nd?\0), 



(4) 



where k is a proportionality constant. The time-dependent 
BCS Hamiltonian Hscsit) is symmetric under the unitary 
fermion parity operator P which transforms fermions if>' (r) 
as P : ij)\(r) — > — ^{r). Thus the initial ground state |0) 



must be an eigenstate of P which has eigenvalues ±1 (since 
P 2 = 1). This is referred to as the ground state having even 
or odd parity. Since P commutes with the BCS Hamiltonian 
HBCs{t) at all times, it must also commute with the uni- 
tary time evolution U. However <v anti-commutes with the 
fermion parity operator P. This leads to a direct contradic- 
tion with Eq.|4] ruling out the possibility s x s 2 = 1. Therefore 
S1S2 = —1, establishing in general the relative — sign for 
non-Abelian statistics. 

Thus there are 2 possibilities for the result of an adiabatic 
exchange 



7i A72 and 72 



-A71 



(5) 



where A = ±1. This operation can be compactly represented 
in the Majorana space as 



(6) 



as claimed. Here the sign in the exponent is determined by 
the path of the adiabatic exchange. The Abelian phase 4> can- 
not be determined by considering only the operator dynamics 
and requires consideration of the ground state wave-function. 
As discussed in the previous subsection, this braid matrix is 
uniquely determined by the Hamiltonian. 



III. T-JUNCTION EXCHANGE 

Consider a junction of three wire segments (A, B, and C) as 
pictured in Fig.[TJ Initially (at time t — 0), a portion of the up- 
per left segment (A) is prepared in the topological regime, re- 
sulting in two Majorana end states labeled at the points a' and 
a. The procedure to exchange these ends through the junction 
takes place in three steps, with step j completed at time t = tj, 
and £3 = t. For a clockwise exchange gates are activated 
(1) to extend the topological superconducting phase into seg- 
ment C; The gates are then deactivated in segment A, result- 
ing in a contraction of the topological region into segment C. 
This procedure is repeated (2) to move the topological region 
into segment B, then (3) back into A. In this process, the two 
ends of the topological region have been exchanged. As noted 
above (Eq.[5]), the resulting transformation must have the form 
{7a'(r) = A 7a (0), 7a (r) = -A 7a ,(0)}, where A = ±1. We 
shall demonstrate this fact for the most general form of the 
junction, and show that the effect of an exchange through the 
junction is the same for Majoranas with the topological and 
with the non-topological phase in the region between them. In 
the process, we shall also show that the unitary transformation 
enacted by an exchange through the junction is independent of 
the overall sign chosen in the definition of any Majorana state. 

In order to determine A, we will examine more closely the 
steps taken in a clockwise exchange. We base this analysis 
on the following three assumptions about the motion of Ma- 
joranas in the wire network: 



1 . All regions of the wire beside the Majoranas have unique 
ground states for both the topological and non-topological 
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that 



B A 



B A 



B A 



FIG. 1: (Color Online) This series of figures shows the interchange 
of two Majorana bound states at the ends of a (solid) topological 
region through the activation and deactivation of the wire segments 
surrounding the junction. Non-topological regions are indicated by 
dashed lines. At the end of step 3, the system has returned to the 
configuration shown at the beginning of step 1, but the endpoints 
have been exchanged. 



phase. That is, aside from the Majoranas at the phase 
boundaries, there are no zero modes in the system. 

2. The process of altering the phase boundaries is carried out 
adiabatically, so that when a topological region is extended 
or contracted from location x to location y, the Majorana 
zero mode operator at its end is transported according to 

7^(2) = £ y xlx(0) for extension 

and 

7 y (t) = C yx ^ x (Q) for contraction. 

3. The process is reversible, implying that £ xy — C yx . 



Without these assumptions, it is impossible to consistently 
perform braiding operations through the junction. For con- 
sistency in expanding within the same segment, we require 



f — f f 

^zx ^zy^yx 



(7) 



so long as either points x and y, or points y and z are in 
the same segment. Likewise, the requirement that Majorana 
fermions remain properly normalized leads to £ xy = ±1 and 

C xy — il. 

Note that Assumption 1 does not preclude the possibility of 
localized fermion bound states, but requires that such states 
are affected deterministically by the passage of a Majorana. 
What is more, it requires that the effect of passing a Majorana 
over such localized states is the same independent of whether 
a topological region is being contracted or expanded. With 
Assumption 2, this is equivalent to the condition that two Ma- 
joranas brought together by the contraction of a single topo- 
logical (or non-topological) region will be the either in the 
ground state or the excited state independent of the location at 
which they are brought together. 

If the parity of a short topological region at x, as determined 
by the fusion channel i^ x _ t ^ x of the two Majoranas at its end 
points, does not match the parity of the ground state when that 
region is eliminated, a quasiparticle will be left behind by the 
elimination of the topological region. This allows us to define 
the local vacuum channel v x for two Majoranas by the form 
of the Hamiltonian 



H = -iT^ t \x,e) lx . e - /x 



(8) 



when they are brought very close together, where r( n ' { ) is the 
coupling constant for two Majoranas with a topological (t) or 
non-topological (n) region between them. We set 



4"' 4) = limsgn (r< n,t >(a?,e)) 



(9) 



where x — e is always further from the junction than x. Con- 
sider the situation in which two Majoranas at points x_ and 
x are endpoints of the same small topological region and are 
transported to points y_ and y closer to the junction. (x_ in- 
dicates a point infinitesimally further from the junction than 
x ) First the right side of the region is expanded from x to y, 
then the left side is contracted from x_ to y_ . Assumptions 1 
and 3 together imply that if the system started in a state such 
that elimination of the topological region would leave the sys- 
tem in an excited state when the elimination happens at x, the 
same must be true if the region is eliminated at y. Likewise, 
if the system began in the ground state of the Majorana inter- 
action Hamiltonian ([8]) at x, it must end in the ground state of 
the corresponding Hamiltonian at y. Therefore, 



C 



•I X- — v x v yC X y 



if points x and y are on the same segment, and 

Cy X _ — ~ v x v yCxy- 



(10) 



(11) 



if points x and y are on different segments. The minus sign 
here occurs due to the fact that j y and j y _ anti-commute, and 
y_ is further from the junction than y. If x and y are endpoints 
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of different topological regions, then the same arguments lead 
to 



C — v n v n C 

<"X-y- — v x v y^yx 

if x and y are on the same segment and 

r - -v n v n C 

^y-X — "i Vy^x-y 

if x and y are on different segments. 



(12) 



(13) 



A. Exchange of Majorana endpoints of a single Topological 
region 

At the first step of the exchange process, the Majoranas are 
moved from segment A into segment C by first expanding the 
topological region to move the Majorana that begins nearest 
the junction (at position a) through it to position c' in segment 
C and then contracting the back end of the topological region 
through the junction to bring the second Majorana from a' to 
c (Fig.[T|. In terms of the Majorana operators, we have 



7c- (h) = £ c ' o 7a(0) 

7c(ii) = C ca '7 a '(0). 



(14) 



At the next step of the process, we use the same procedure to 
move from segment C to segment B, leading to 

Jb'(h) = £b'clc{h) 

7b(*a) =C 6 c'7c'(*i)- (15) 
Finally, the Majoranas are moved back into segment A, with 

la'{h) = ^a'blbih) 

la{h)=C ab , lb >{t 2 ). (16) 

In total, we have that 

7a' (t) = £a'(Ac'£c'a7a(0) 

7 a (r) = C ab >£ b > c C ca >7a'(0)- (17) 
Using Eq. Q, we can reduce this to 

7a' (t) = £a'a-£a-bCbc£ca7a{0) 

= C a _ o ^*x*7a(0) 

7o( r ) = Cab£bcC C a-C a _ a ,-f a ,(0) 

= -Ca_a'^X*7a'(0), (18) 

where we have defined 

X* = v l a £ a _ h C hc £ ca (19) 

and used Eq. ( fTI) to permute the indices in the second equa- 
tion. Note that is defined uniquely for the junction, since 

V a £a-bCbc£ca — ^ b £b-c(-'ca£ab — V c £ c _ a C a b£bci 



which can be shown using Assumption 2 and Eq. (Hi. Fur- 
thermore, due to Eq. ^ x* is independent of the locations 
of points a, b and c. It is only important to the definition of 
X* that each of these three points be on a different one of the 
segments connected by the junction. 
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FIG. 2: (Color Online) This series of figures shows the interchange of 
two Majorana bound states at the ends of a (dashed) non-topological 
region through the activation and deactivation of the wire segments 
surrounding the junction. 



B. Exchange of Majorana endpoints of separate Topological 
regions 

Suppose now that the two Majorana bound states that we 
wish to exchange begin as endpoints of different topological 
regions, as shown in Fig. [2] We begin the exchange process by 
contracting the topological region from segment A into seg- 
ment C, so that one of the Majorana bound states moves from 
point a to point c. We then move the remaining Majorana 
across the junction by expanding the other topological region 
from point a' to point b in segment B. This leads to 



7c (ti) 
lb{h) 



Ccala{0) 
£ba'la'{0)- 



(20) 



The time ti here is the time at which step i in Fig. [2] is com- 
pleted. We are now faced with something of a dilemma: in or- 
der to move the Majorana from segment C back into segment 
A, we must first extend the topological region from point c 
to the junction, and then contract it into segment A. During 
this process, we bring three topological regions together at the 
junction, a step that was unnecessary in the previous type of 
exchange. 

Due to the ambiguity in which of the three segments the 
junction point lies on, we cannot use our previously estab- 
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lished C and £ moves to carry out this process. Instead, we 
avoid the ambiguity by introducing a new factor U to describe 
this type of move (expand into junction and contract away). 
Then 



7o'(*3) = Ua'cldh) 



(21) 



The requirement of proper normalization for 7 a and ~f c still 
holds, so U a ' c — ±1. Reversibility of the process leads to 
U a i c = U ca i . 

Continuing with the exchange, we contract the topological 
region from segment B into segment C, then extend it back 
into segment A to bring the remaining Majorana to point a. 
That is, 



In total, we have 



7a' 0) 

7a (r) 



t^a' c^ca 



7a (0) 
£acC c b£ba'Ja<{0)- 



(22) 



(23) 



By the argument of Sec. II B and given Assumption 1, the two 
prefactors above must differ by a — sign. That is, it must be 
that 

U a 'c = —C c b£ba' (24) 

Because the points a', b and c are arbitrary, this equation 
must hold whenever all three points are in different segments 
around the junction. 

As we have above, we now separate the effect of the ex- 
change process two factors, one describing motion along A 
and one a property of the junction itself. Defining x™ by 

X n = v%£ ab C bc £ ca _ ) (25) 

we have that 



7a' (t) 
7a (t) 



£a_a'<X"7a 
^a_a' u a A /a' ) 



(26) 



where we have used Eq. ( 13 i to permute indices where nec- 



essary. As in the case of Majoranas connected by a topolog- 
ical region, the value of x is consistent when under transfor- 
mations that rotate the three segments (A — >• B — > C), and 
antisymmetric under the exchange of any two segments. As 
with x*, X™ i s independent of the points chosen on the three 
segments for the representation given by Eq. ( |25] l. Because 
of these properties, we shall refer to \ as the chirality of the 
junction, and show that x = X™ = X*- 

In order to see this, it is instructive to consider another pro- 
cess, in which two Majoranas connected by a non-topological 
region are transformed into two Majoranas connected by a 
topological region (Fig. [3]). It is shown by this process, as 
well as by the arguments leading to Eqs. (Hi and ( 13 1 that 



When combined with Eq. (24]i, this gives 

Cbc£cci- 



v a l v t b £ b _ c C c 



U ba _ = -v a l v t b e b _ c c ca . (27) 
ives 

(28) 

Using this with Eq. (20i and Eq. (25 1 shows that \ n = x*- 
That is, the effective junction chirality is independent of 
whether the Majoranas involved in the exchange are endpoints 
of the same or different topological regions. 
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FIG. 3: (Color Online) Majoranas at the endpoints of different topo- 
logical regions are transformed into Majoranas at the endpoints of 
the same topological region by a three step process. The third step 
requires the use of the factor Ub a _ as in Eq. 424b. 



C. Braid transformation 

We can combine the two cases for Majorana braiding into 
a single notation, writing the unitary time evolution operator 
for the entire process as 



U x {r) 



(29) 



where M. is the sign (either C or £ ) required to move the Ma- 
jorana at a' to the point a_ close to a. v a is the vacuum chan- 
nel as determined by the local Hamiltonian of the form |[8). 

This operator has the advantage of invariance under the 
sign ambiguity inherent in the Majorana fermion description. 
While the time-evolution of a Majorana fermion leads to a 
definite result 7(7-) = UUt)^(Q)U x {t) for, it is important 
to remember the following caveat about the definition of the 
Majorana operators themselves. If j(t) is a zero energy op- 
erator (in the sense that [HBcs(t),lj(t)] = 0), then so is 
7j(<) = vy(t) for v = ±1. This is essentially a phase am- 
biguity in the definition of a Majorana fermion (note that a 
general phase e %e is not allowed since 7' (t) = 7(t)). How- 
ever, in redefining any particular Majorana operator involved 
in the exchange, one must also change the sign of any vacuum 
channel, contraction factor or expansion factor in which that 
Majorana is involved. For example, in redefining 7 a — > —-f a , 

one must also change — > —Va' n \ C xa — s- —C xa and 
£ xa — > —£ xa in order to maintain consistency with the defi- 
nition of other Majorana operators. Because the definition of 
X n (or %*) contains factors involving each Majorana operator 
twice, x is invariant under this transformation. Likewise, if 
one changed j' a — > — 7^, U x as defined above remains invari- 
ant because A4 a -a' would also change sign. 
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IV. DISCUSSION 



A. Comparison of Junction Chiralities 



It is important to note that the chirality of a given junction 
is not set a priori by the general analysis considered here. In- 
stead, it is set by the underlying microscopic parameters and 
the exact method by which Majoranas are moved through the 
junction. In particular, a junction's chirality may be altered 
by a local defect that binds a fermion in only one of the two 
phases. Such a defect causes the Majorana bound states to 
emit or absorb a fermion whenever they pass by the defect, 
changing the sign of each Majorana. If the defect is located 
at the junction, it will be passed three times by the Majoranas 
during an exchange. Each Majorana would therefore acquire 
an additional minus sign during the exchange, changing the 
effective chirality of the junctionP^l 

In any case, different junctions within a wire network do not 
necessarily have the same chirality. To understand the effect 
that this might have on a computational algorithm, consider 
the situation in which a pair of Majorana fermions 72 and 73 
are exchanged clockwise via first one junction with x — Xi 
then through a second junction with x — X2- If Xi 7^ X2, one 
of these exchanges is effectively a counter-clockwise, rather 
than a clockwise, braid. In this case, the net effect is no braid 
at all, rather than a double braid. 

More explicitly, we shall examine the effect of this manip- 
ulation when 72 forms a qubit with ji, and 73 forms a qubit 
with 74. Then the fusion channel of 71 and 72 evolves as 



' ba 



(30) 



That is, the 71, 72 qubit flips do to this double exchange only 
if the chiralities of the two junctions being used are the same. 
Likewise, the 73, 74 qubit also flips only if xi = Xi- Along 
with a method of measuring the qubit states (such as the frac- 
tional Josephson effect as propose in Ref. 1), this double- 
exchange test may be used to determine the relative chirality 
of different junctions in a wire network. 



B. Multiply connected networks 

It is worth noting that while the result of a braiding of two 
Majoranas is uniquely determined by the properties of the un- 
derlying wire network, the braid result may be dependent on 
the path taken through the wire system. Consider, for instance, 
the possibility that segments B and C in the T-junction are 
connected via a loop of wire. It would then be possible to in- 
terchange the Majoranas at positions a and a 1 in Fig.[T]clock- 
wise by taking them from segment A to segment B then trans- 
porting them around the loop to segment C, and finally back 
through the junction into A. The result of this transportation 



7a'(r) = £a'c£cb'£b'aja(0) = C a _ a> V^X^-l a{fy 
7a (t) = Cac'CdbCbaHa'i®) = -C a _ a > V* X^7a (0) ,(3 1) 





£ cb £ ba 



U ac £ cb £ ba 



FIG. 4: (Color Online) This series of figures shows the transport 
of a Majorana fermion around a loop in a multiply connected wire 
network. Each figure indicates the total phase factor acquired in the 
motion of the Majorana to that point. Since U a c = ~C c b£ba the total 
phase factor acquired by the Majorana is C — —£ c b£bc- 



where L = —£ c b£bc = —C c bCb c , and where the £ c b> indicates 
expansion around the loop from b to c, rather than through 
the junction. That £ c b£bc = C c bCb c may be shown using As- 
sumption 2, Eq. |7]) and Eq. (111. The result here differs from 
Eq. ( 18 1 by the 'loop factor' C. This factor is independent of 
the points b and c used in the definition, as well as being even 
under the exchange of segments B and C. C = — 1 may be 
thought of as indicating the existence of a flux (real or effec- 
tive) through the loop in the network. 



Consider Majorana fermions at points a' and a near a loop 
in the wire network as pictured in Fig. [4] Transporting the 
Majorana around the loop and back to its original position 
yields a transformation 



7a ( 



£7a(0). 



(32) 



Since 7 a < is unaffected by this process, if C = — 1 the fermion 
parity operator in the (a', a) (proportional to t~f a >j a ) changes 
sign. This indicates that if there were an odd number of 
fermions in the (a', a) region before the process began, then 
there will be an even number afterward and vice versa. This 



may appear to contradict the conclusions of Sec. II A how- 
ever, the conclusion that the fermion parity is conserved in 
the motion of Majorana fermions was dependent upon the as- 
sumption that Hbcs( t ) = Hbcs(Q)> which is not the case 
here. In fact, the loop has switched phase. The extra fermion 
parity may therefore be found in the loop, either in a localized 
bound state or in the alteration of the ground state of the loop 
from even to odd parity if a real external flux is present. In 
this sense, C is seen to indicate the presence or absence of an 
effective flux through the loop. 
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C. Applicability 

The results of this paper apply to any one-dimensional wire 
system supporting Majorana fermions at the boundaries be- 
tween two phases for the system, so long as the assump- 
tions laid out in Sec. [In] hold true. The details of the un- 
derlying implementation of such a system, e.g. in semicon- 
ductor/superconductor heterostructures with large spin-orbit 
coupling 1 , will determine the values of the factors \ and L 
that we have identified here for junctions and loops, respec- 
tively. Our results establish a concrete framework for describ- 
ing and tracking the non-Abelian transformations made by 



moving Majorana fermions through a quantum wire network. 
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